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Abstract 

A model was proposed to account for a new kind of history effect in the 
transport measurement of a sample with inhomogeneous flux pinning coupled 
with flux creep. The inhomogeneity of flux pinning was described in terms of 
alternating weak pinning (lower jc) and strong pinning region (higher jc). The 
flux creep was characterized by logarithmic barrier. Based on this model, we 
numerically observed the same clockwise V-I loops as reported in references. 
Moreover, we predicted behaviors of the V-I loop at different sweeping rates of 
applied current dl/dt and magnetic fields B^, etc. Electric transport measure- 
ment was performed in Ag-sheathed Bi2_xPb2:Sr2Ca2Cu30jy tapes immersed 
in liquid nitrogen with and without magnetic fields. V-I loop at certain dl/dt 
and Bq was observed. It is found that the area of the loop is more sensitive 
to dl/dt than to Bq, which is in agreement well with our numerical results. 
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I. INTRODUCTION 



Recently the so-called history effect (HE) of vortex matter has been frequently observed 
by a variety of methods such as electric transport measurement, dc magnetization hysteresis, 
ac susceptibility and torque techniques, and drawn much more attention |ll|Pj3|J^,p|,|6|,[7|,p|,|9|,p!0 



HE refers to that critical current density jc of a superconductor has different values at a same 
field (or temperature) for either a thermal, magnetic or current cycles. The earlier observed 
hysteresis of R-T curves in the thermal cycle in fact is a kind of HE caused by first order 
phase transition of vortex matter (melting or freezing) [ IT | . Different mechanisms have been 
proposed to account for HE ]l|,|^J^J|,|,^^||,||,|lO| . And a typical explanation contains two 
assumptions: relatively weak flux pinning and uniform small quenched disorders. Then HE 
is explained with a theory of disorder-order transition of vortex matter. The disordered 
vortex matter pinned more strongly (higher jc) is supercooled to low temperature where 
it is in metastable state and then a process such as a bias current anneals the metastable 
disordered state into an ordered and stable one with lower jc. An important point is that 
flux creep in the models is omitted for the low temperatures. In this case, the property of 
vortex matter is governed only by the competition between interactions of vortices- vortices 
and vortices-quenched disorders. 

However, a different kind of HE has been observed [|T2|. The hysteresis loop of V-I 
curve of polycrystalline Bi2-xPbxSr2Ca2Cu30y (Bi2223) at liquid nitrogen temperature was 
measured there. See also Fig. 6 in this paper. It is clear that such a V-I loop is a new kind 
of HE and very different from the V-I loop in low temperature superconductors reported 
before The difference between these two kinds of HE is that the direction of the V-I 
loop is anti-clockwise in reference 1 while clockwise in reference 12 and this paper. Despite 
the anti-clockwise loop has been explained very well, the clockwise one has not yet been 
understood up to our knowledge and interested us very much. 

Note that this new kind of HE takes place in a system different from those with uniform 
and weak quenched disorders at low temperatures. First, at higher temperatures, such as 
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liquid nitrogen temperatures, thermal fluctuation should be important and flux creep could 
influence the property of vortex matter apparently [0 . In fact the influence of flux creep on 
the V-I curve and transport property of Bi2223 have been conflrmed by earlier experimental 
and numerical results p!^JT5|JT6|jr7|JT8|] . Second, the quenched disorders or flux pinnings are 
inhomogeneous in polycrystalline high temperature superconductors (HTS). There are weak 
links and grains in sintered samples such as polycrystalline Bi2223. And the flux pinning 
strength of the weak links at grain boundaries is very different from that of the grains [jl9| . 
In view of these two points, any models concerning the system with homogenous quenched 
disorders without flux creep could not be used to understand this novel hysteresis behavior 
in V-I curve. 

A theoretical model describing irreversible electric-magnetic behaviors of a system with 
inhomogeneous quenched disorders has already been proposed |2^ . Unfortunately, flux creep 
was omitted there. Moreover, this model has not been used to explain any HE, including 
the clockwise V-I loop. Experimentally, A.M. Campbell's group did many important works 
on inhomogeneous superconductors. They did observe the hysteresis V-I curves before. 
However, only a qualitative explanation was given. And what interests us now is how to 
simulated the clockwise V-I loop. 

Considering flux creep and inhomogeneous flux pinning, we proposed a new model to 
explain this new kind of HE. We carried out numerical simulation with this model. Upon ex- 
plaining the clockwise loop of V-I curve reported in reference 12, we further predicted some 
new phenomena. To test the predictions, electric transport measurement was conducted 
for characteristic V-I curve of silver-sheathed Bi2_a;Pba;Sr2Ca2Cu30j; tapes (Ag-Bi2223) im- 
mersed in liquid nitrogen with or without applied flelds B^. Clockwise loop of V-I curve 
at certain dl/dt and B^ was observed and is dependent on dl/dt and B^, conflrming the 
numerical result. 
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II. MODEL 



The model consists of two basic assumptions: there exists flux creep in a sample and 
the flux pinning is inhomogeneous. The complex inhomogeneity of flux pinning is simplified 
by periodic strong pinning region (S) with large jc (jcs) and weak pinning region (W) with 
small jc iicw), see Fig.l. The flux diffusion is characterized by collective creep of vortex glass, 
which has been discussed widely To simplify the calculation, we use logarithmic U(j) 
relationship, which is the special case of vortex glass model with very small glass exponent 
/i 

We should emphasize that our model is not a simplified version of the "brickwall" model. 



which has been widely discussed before ||2^. There are some resemblances between them. 
However, there are at least two points that distinguish our model from the previous one. 
First, the "brickwall" model is a static critical state model without flux creep. Obviously, in 
any static models, the transport properties studied here will be independent of the sweeping 
rate of applied current dl/dt, contrasting with the experimental observation. Second, in the 
"brickwall" model the S and W regions are connected in series along the direction of current. 
While in our model, they are connected in parallel. 

In our one- dimensional model, flux lines firstly enter the system from the surfaces into 
the weak pinning channel (W) and then diffuse into the strong pinning region (S). So W 
surrounds S. In other words, S is inside the slabhke sample. 

As an example, the model may be proper to polycrystalline samples of HTS, such as 
silver-sheathed Bi2_a;Pba;Sr2Ca2Cu30j, (Ag-Bi2223). The grain boundaries in polycrystalline 
HTS form weak link network of flux pinning and thus can be considered as W, whereas the 
intra-grain region is S. It is well known that the weak link network suppresses substantially 
the critical current density of HTS such as Ag-Bi2223. When current or magnetic field is 
applied, flux density will enter inter-grain (W) at first and then penetrate gradually into the 
intra-grain (S) by means of flux creep. Another example the model may be applicable to is 
the low temperature superconductors, such as NbsSn metallic compounds, where flux lines 
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may be pinned by grain boundary (S) whereas the intra-grain is weaker pinning region (W) 

Certainly, one cannot expect that numerical result by such a simplified phenomenological 
model can quantitatively represent experimental data of a sample. For example, the one- 
dimensional geometry assumption in the model is not fulfilled for a real sample. And, the fiux 
pinning is more complex by far than the two model parameters, ics and icw- Nevertheless, 
our model is a new one accounting for the experimental clockwise V-I loop, which was 
observed very recently. The effectiveness of this model was demonstrated by the electric 
transport measurement of hysteresis V-I loop. 

III. SIMULATION 
A. Basic Equation 

Consider a slab with infinite length along y-axis, thickness d along the x-axis and width 
w along the z-axis. The current is applied along the y-axis. In view of w>>d, we focus 
on the one- dimensional case for simplicity. The non-linear diffusion equation describing 
the macroscopic B or j is well known p5| , pB| . With a logarithmic dependence of U on j: 
U{j) = Uq In \jc/j\, the diffusion equation of flux can be written as: 

dB vq d 

'dt ~ (/ioJc)"+' 

where n = U^/kT and Vq = uuOm- Vq is the attempt velocity of the thermal-activated 
vortex motion, u is the hopping distance and Um is the microscopic attempt frequency. 

B. Boundary and initial conditions 

In the V-I curve measurement, the current is always applied with a certain sweeping 
rate dl/dt. The boundary condition of Eq.l can be obtained from the current conservation 
equation: ^ Jq jdxdz = ^. Substituting dB/dx = —fioj into it, one comes B{d,t) — 
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dB 
dx 




(1) 



B{0,t) = —fiol/w. Due to antisymmetry, one has B{d,t) = —B{0,t), thus the boundary 
condition is 

B{0,t) = ^=(l^^]t (2) 

^ ' ' 2w \2wdt ^ ^ 



As for the initial condition, there is no current at t=0 in the sample, hence 

B{x,0)=Q (3) 

C. Numerical method and the choice of parameters 

Such kind of non-linear diffusion equation can be numerically solved with the finite 
difference method. First, the temporal and spatial variances are discretized with different 
steps respectively. Since the temporal and spatial steps are very important for the stability 
of the calculation, we take pains to choose suitable steps to avoid a divergence problem. 
Second, the equation is discretized and rewritten in an implicit difference scheme, which is 
better than the explicit difference scheme in the stability of calculation. Finally, based on the 
initial condition, the calculation by iteration will obtain the temporal and spatial magnitude 
of physical quantity, namely B(x,t), j(x,t) and E(x,t). The voltage can be obtained by 
integration of E(x,t). Note that the boundary condition plays a role of constraint in every 
time step during current ascending and descending, see Eq.2. 

It is easy to see that U{j) = Uolia\jc/j\ combined with E = Eoexp{—U/kT) leads to 
the power dependence of the E-j curve E = £'o(j/jc)" with n = U^/kT. This power law E-j 
curves are observed frequently in transport measurements and can be used to determine the 
parameter n [|T5|. Certainly, n is dependent on flux pinning strength, magnetic field and 



temperature |T8| , |27[] . The weaker flux pinning and higher temperature bring about smaller 
n. In fact, many experiments including the present work show that for the Ag-Bi2223 in 
liquid nitrogen temperature, n = 6 is a typical value. 
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As for the velocity of the thermal- activated vortex motion and the critical current den- 
sity, typical and reasonable magnitudes were employed based on the earlier theoretical and 
experimental work, for example we took vq = uujm = ^rn/s ( if u ~10~^m, ujm~ and 
jcw — 2 X lO^A/m^. It should be pointed out that our numerical result is not sensitive to 
the choice of n, vq and jcW- 

D. Numerical results and discussions 

1. Effect of jcs/jcW ratio on the hysteresis loop 

Since jc represents the pinning strength, the ratio jcs/icw can be considered as a param- 
eter reflecting the ratio of the two pinning strengths. Shown in Fig. 2 is our numerical result 
of the dependence of HE on ics/icW at fixed icw, ds/dw and dl/dt. It is very clear that the 
V-I loop is clockwise as reported in reference 12. To understand how jcs/jcW affects the V-I 
loop, the corresponding current distributions in sample were also calculated. 

Note that the velocity of flux diffusion, namely the velocity of current diffusion in trans- 
port measurements, can be written as: v — voexp{—U/kT) oc (j/jc)"' for logarithmic U(j). 
Hence, for the same j, the larger the jc, the smaller the v. Generally speaking, the average 
speed of current diffusion in S region ( ) is smaller than that of W region {u^). In our 
calculation, jcW is fixed, i.e. is fixed. By changing jcs/hw ratio, namely jc5, we adjust % 
only. 

For a large ratio {ics/icW—^O), i-e. high inhomogeneity of flux pinning, is so small 
comparing with u^. And in such a case, there is no enough time for current diffusion in 
response to both the ascending branch of V-I curve (Fig.2.1(b)) and the descending one 
(Fig. 2. 1(c)). Consequently, no obvious V-I loop can be detected (Fig. 2. 1(a)). 

For a medium ratio (jcs/hw—^:), i.e. a medium inhomogeneity of flux pinning and a 
larger it is shown that the area of the loop is very large (Fig.2.2(a)). In this case, 
although current has already penetrated into S region during the ascending branch (Fig. 2. 2 
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(b)), much more current penetrates into S during the descending branch (Fig.2.2(c)). That 
is to say, the current diffusion of S region has more time in response to the descending branch 
than to the ascending one. As a result, the V-I loop is obvious. 

For a small ]cs/}cw ratio, low inhomogeneity of flux pinning, vl is very large and ap- 
proaches u^. And there is little difference between the distributions of current during the 
two branches (Fig. 2. 3(b), (c)). The minor difference is only seen when the current is small 
and the corresponding voltage is too small to be detected. Hence, there is almost no hys- 
teresis loop (Fig. 2. 3 (a)). Especially, at the uniform case {]cs/']cW =1, ^ = there will 
be no loop at all. 

2. Effect of dg/dw ratio on the hysteresis loop 

With a fixed total thickness, the spatial distribution of S and W region, i.e. the thickness 
ratio of them is also important for the V-I loop. The numerical V-I loops with different 
ds/dvK at fixed dl/dt and ]cs/']cw were displayed in Fig.3. It is found that the smaller the 
d^/dvK, the smaller the area of the loop. 

For d^/diy =0.857 and other parameters given here, the thickness of S region i.e. the 
distance for current diffusion in S region is so long that current diffusion in S region has 
less time to respond the ascending branch (Fig. 3. 1(b)) than to respond the descending one 
(Fig. 3. 1(c)). So the area of the V-I loop is large (Fig.3. 1(a)). 

For d^/dvK =0.222, S region is thinner, i.e. the distance of current diffusion in S region 
is shorter, resulting that comparing with the <ls/^w ==0.857 case, the current diffusion in 
S region has more time in response to the ascending current (Fig.3. 2(b)). In other words, 
for such parameters given here, the response time of current diffusion in S region to the two 
branches of the V-I loop shows httle difference. As a result, the area of the V-I loop becomes 
smaller (Fig. 3. 2 (a)). 

As for ds/<^w =0.105, S region is too thin and the distance of current diffusion in S 
region is too short. Hence there is enough time for current diffusion in response to both 
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current ascending and descending at the given parameters (Fig. 3. 2(b), (c)). That is to say, 
there is no obvious V-I loop (Fig. 3. 3(a)). As ds/dw approaches 0, namely the uniform flux 
pinning case, the hysteresis loop will disappear thoroughly. 

3. Current sweeping rate (dl/dt) dependent V-I loop 

The numerical result of current sweeping rate dependent V-I loops at fixed ]cs/]cW and 
ds/dvF was illustrated in Fig.4. 

It is seen that if dl/dt is very small, such as O.OlA/s, fiux lines will have enough time to 
respond the changing applied current and thus can penetrate into both W and S regions. The 
difference of current distribution during the ascending (Fig. 4. 1(b)) and descending branch 
(Fig.4. 1(c)) is very small. As a result, the area of the V-I loop is also small (Fig.4. 1(a)). 
As dl/dt reaches a moderate value, say dI/dt=2A/s, the difference between the current 
distributions during the ascending and descending branch is larger. And the area of the V-I 
loop becomes larger too(Fig.4.2(a)). Nevertheless, when dl/dt is too large, say dI/dt=10A/s, 
the area of the V-I loop decreases dramatically (Fig. 4. 3 (a)). It is easy to understand that in 
such a case, comparing with dl/dt, vl is so small that the current cannot penetrate into the 
S region during both the ascending and descending branch. And the corresponding current 
distributions during the two branches do not show visible difference (Fig.4.3(b), (c)). So 
there is no obvious V-I loop. 

Now, it is clear what the effect of dl/dt on the V-I loop means. With fixed distance and 
average speed of the current diffusion, to change dl/dt is equivalent to shift our observing 
time window. When we choose time window by using dl/dt as lOA/s, the observing time 
is so early that there is no enough current inside S region by means of diffusion even the 
applied current sweeping has been finished. Of course there is no detectable V-I loop. If we 
use dI/dt=0.01A/s, we shift our time window so late that there are always enough current 
diffused into S region during both the ascending and descending branch. It is natural that 
no obvious V-I loop exists at all. Only when we choose a proper observing window, i.e. a 
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proper dl/dt such as 2A/s, is there much more current in S region during the descending 
branch than the ascending one. In other words, only in such case, can the current diffusion 
in S region has more time to respond the descending branch than the ascending one. And 
the V-I loop will be obvious. 

In conclusion, if there is more time for current diffusion of S region to respond the 
descending branch than to respond the ascending one, hysteresis loop in V-I curve, i.e. the 
HE will be obvious. On the other hand, if the current diffusion time of S region is similar 
or the same for the ascending and descending branch, there will be no obvious V-I loop, 
namely no HE. The observation of HE can be controlled by adjusting dl/dt ( the observing 
time window ) or ds/'^w ( the distance of current diffusion ) or ]cs/]cw ( the average speed 
of current diffusion ). 



4- Magnetic field dependent V-I loop 

It is well known that jc is dependent on the magnetic field Bq. Commonly, jc decreases 
with increasing B^. Therefore, to study the effect of jc on the V-I loop is equivalent to study 
the effect of B^ on it. For the inhomogeneous pinning case, ]cw is more sensitive to B^ than 



jcs, especially when B^ is not very high ||T9|. Thus, we simulated the V-I loops in different 
applied fields B^, namely different ]cw at fixed jcs, dl/dt and dg/diy- Shown in Fig.5 are 
the numerical results. It is clearly seen that the loop shifts towards smaller current with 
decreasing ]cw while the area of the loop is insensitive to ]cw- In other words, the area of 
the loop is weakly affected by magnetic field. 

As mentioned above, the existence of HE depends on that the current diffusion of S region 
has different responding time for the two sweeping branches of V-I loop. Only changing ]cw, 
i.e. Bq, would not affect the responding time of S region. As a result, the area of V-I loop 
will not change acutely. 
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IV. EXPERIMENTAL 



Because our numerical simulation has obtained more results than the experimental one 
reported in reference 12, such as the effect of dl/dt on the area of loops, we measured V-I 
curves of Ag-Bi2223 samples in various applied magnetic fields Bo with different dl/dt to 
confirm the numerical prediction. It has been pointed out above that the Ag-Bi2223 sample 
may be proper to test our model. 

The Ag-Bi2223 samples cut from long silver sheathed tapes come from the National 
Center for R&D on Superconductivity of China. All the samples are the c-axis textured. 
The c-axis is perpendicular to the wider surfaces of the samples. Including the outer silver 
sheath, the average size of the samples is approximately: 4mm x 1mm x 4cm. Between each 
two neighbor points of four (two voltage contacts and two current contacts) is 1cm in length. 
The transport measurement was carried out for the sample immersed in liquid nitrogen. The 
external magnetic field with range ~ 200 Gs was applied in the superconducting state. The 
current range of the supply source is ~ 40 A and sweeping rate range 0.02 ~ 1 A/s. The 
measurements were made in the case of Ba parallel to the c-axis of the samples. More details 
of the sample preparation and the measurements can be found in reference 18. 

The V-I loops were measured by ascending current first and then descending it. Displayed 
in Fig. 6 are the experimental V-I loops with different dl/dt and fixed B^. The spottiness 
of the experimental data came from the limited picking rate of our measurement system. 
As reported in reference 12, clockwise loop of V-I curve was observed. It is clearly seen 
that dl/dt affected not only the appearance of V-I curve but also the area of the loop. For 
example, the area of loop with dl/dt = 1 A/s is larger than the area of loop with dl/dt — 
0.5 A/s, which is in agreement well with the numerical result. 

Shown in Fig. 7 are the experimental V-I loops in different applied fields Ba at a fixed 
dl/dt. It is found that the area of the loop is insensitive to magnetic field Ba whereas 
the transport critical current density jc decreases with increasing applied field very sensi- 
tively. This experimental dependence confirms the above numerical results and thus supports 
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strongly our model. Furthermore, the agreement between the numerical and experimental 
V-I shifting by applied field implies that the critical current density jc of a sample with 
spatially non-uniform flux pinning strength, such as the silver sheathed Bi2223 tapes, is 
governed mainly by the weaker flux pinning region (the region with smaller jc). Hence, for 
these HTS materials, which have potential for technical applications, the point to enhance 
their critical current density is to improve the flux pinning strength of weak pinning region. 

We are unable to measure more V-I curves to test all the numerical curves for the limited 
ability of our instrument. For example, dl/dt cannot be too small because Ohm heat at the 
current contacts will destroy the sample. 

V. SUMMARY 

We have proposed a model to account for a new kind of history effect in an inhomogeneous 
flux pinning superconductor at high temperatures. The model consists of alternating weak 
and strong flux pinning regions whose strength was depicted by different critical current 
densities jcW and icS, respectively. Based on the model, our numerical simulation successfully 
observed a hysteresis loop in the characteristic V-I curve as reported in references. Our 
numerical results also predicted some new characteristics of the V-I loop at different sweeping 
rates of applied current dl/dt, applied magnetic flelds B^, the ratio of the two regions' pinning 
strength and their thickness. 

Physically, due to the difference of flux pinning strength in S and W regions, the average 
speeds of flux (or current) diffusion in the two regions are different too. If dl/dt (or ds/dw, 
ics/icw) is proper, the V-I loop, i.e. the HE is obvious. Otherwise, no obvious V-I loop can 
be observed. 

To conflrm the numerical results, electric transport measurement was conducted to get 
the characteristic V-I curve of Ag-sheathed Bi2-a;Pba;Sr2Ca2Cu30j/ tapes immersed in liquid 
nitrogen with or without applied flelds. Clockwise V-I loop was observed to be dependent 
on dl/dt and B^. Hence, we presented a new possible mechanism accounting for the new 
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kind of HE in an inhomogeneous flux pinning superconductor. 
Acknowledgment 

The work is supported by the Ministry of Science and Technology of China (NKBRSF- 
G1999-0646) and NNSFC (No.l9994016). 



13 



REFERENCES 



[1] W. Henderson, E. Y. Andrei, M. J. Higgins and S. Bhattacharya, Phys. Rev. Lett. 77, 
2077 (1996); Z. L. Xiao, E. Y. Andrei, P. Shuk and M. Greenblatt, Phys. Rev. Lett. 85, 
3265 (2000). 

[2] Y. Paltiel, E. Zeldov, Y. N. Myasoedov, H. Shtrikman, S. Bhattacharya, M. J. Higgins, Z. 
L. Xiao, E. Y. Andrei, P. L. Gammel, D. J. Bishop, Nature (London), 403, 398 (2000). 

[3] S. Kokkaharis, P. A. J. de Groot, S. N. Gordeev, A. A. Zhukov, R. Gagnon and L. Taillefer, 
Phys. Rev. Lett. 82, 5116 (1999). 

[4] S. Bhattacharya and M. J. Higgins, Phys. Rev. B 52, 64 (1995). 

[5] S. O. Valenzuela and V. Bekeris, Phys. Rev. Lett. 84, 4200 (2000). 

[6] S. S. Banerjee, N. G. Patil, S. Saha, S. Ramakrishnan, A. K. Grover, S. Bhattacharya, 
G. Ravikumar, P. K. Mishra, T. V. Chandrasekhar Rao, V. C. Sahni, M. J. Higgins, E. 
Yamamoto, Y. Haga, M. Hedo, Y. Inada, and Y. Onuki, Phys. Rev. B 58, 995 (1999). 

[7] R. Wordenweber, P. H. Kes, and C. C. Tsuei, Phys. Rev. B 33, 3172 (1986). 

[8] S. S. Banerjee, N. G. Patil, S. Ramakrishnan, A. K. Grover, S. Bhattacharya, P. K. 
Mishra, G. Ravikumar, T. V. Chandrasekhar Rao, V. C. Sahni, M. J. Higgins, C. V. 
Tomy, G. Balakrishnan, and D. Mck. Paul, Phys. Rev. B 59, 6043 (1999). 

[9] S. S. Banerjee, S. Ramakrishnan, A. K. Grover, G. Ravikumar, P. K. Mishra, V. C. Sahni, 
C. V. Tomy, G. Balakrishnan, D. Mck. Paul, P. L. Gammel, D. J. Bishop, E. Bucher, M. 
J. Higgins, and S. Bhattacharya, Phys. Rev. B 62, 11838 (2000). 

[10] A. A. Zhukov, S. Kokkaharis, P. A. J. de Groot, M. J. Higgins, S. Bhattacharya, R. 
Gagnon, and L. Taillefer, Phys. Rev. B 61, R886 (2000). 

[11] H. Safar, P. L. Gammel, D. A. Huse, D. J. Bishop, J. P. Rice and D. M. Gingsberg, Phys. 
Rev. Lett. 69, 824 (1992). 

14 



[12] P. Usak, L. Jansak and M. Polak, Physica C 350, 139 (2001). 

[13] A. I. Larkin and Y. N. Ovchinnikov, J. Low Temp. Phys. 34, 409 (1979). 

[14] L. P. Ma, C. H. Li, R. L. Wang and L. Li, Physica C 279, 143 (1997). 

[15] S. Y. Ding, C. Ren, X. X. Yao, Y. Sun and H. Zhang, Cryogenics 38, 809 (1998). 

[16] P. Zhang, C. Ren, S. Y. Ding, Q. Ding, F. Y. Lin, Y. H. Zhang, H. Luo and X. X. Yao, 
Supercond. Sci. TechnoL 12, 571 (1999). 

[17] S. Y. Ding, H. Luo, Y. H. Zhang, Q. Ding, F. Y. Ling, P. Zhang, X. F. Wu, L. Qiu and 
X. X. Yao, J. Superconductivity 40, 29 (2000). 

[18] Y. H. Zhang, H. Luo, X.F. Wu and S. Y. Ding, Supercond. Sci. TechnoL 14, 346 (1999). 

[19] Ph. Vanderbemden, R. Cloots, M. Ausloos, R. A. Doyle, A. D. Bradley, W. Lo, D. A. 
Cardwell, A. M. Campbell, IEEE Trans. Appl. Supercond. 9, 2308 (1999). 

[20] R. Wordenweber, Rep. Prog. Phys. 62, 187 (1999). 

[21] D. S. Fisher, M. P. A. Fisher, and D. A. Huse, Phys. Rev. B 43, 130 (1991). 

[22] M. V. Feigel'man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Phys. Rev. Lett. 
63, 2303 (1989). 

[23] J. Mannhart, C. C. Tsuei, Z. Phys. B 77, 53 (1989). 

[24] H. UUmaier, Irreversible Properties of Type II Superconductors, Springer- Verlag 89-94 
(1975). 

[25] G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. 
Mod. Phys. 66, 1125 (1994). 

[26] V. M. Vinokur, M. V. Feigel'man and V. B. Geshkenbein, Phys. Rev. Lett. 67, 915 (1991). 

[27] Y. H. Zhang, Z. H. Wang, H. Luo, X.F. Wu, H. M. Luo and S. Y. Ding, J. Phys. Condens. 



15 



Matter 13, 2583 (2001). 
Figure Captions: 

Fig.l. A schematic sketch of the one-dimensional model: a superconducting slab consisting 
of periodic strong (S) and weak (W) pinning regions. Vortex lines enter the system firstly 
from the surfaces into W region, and then diffuse from the W into S regions. So the W 
region surrounds the S region. 

Fig. 2. Numerical results, showing the effect of ]cs/]cW ratio on the V-1 loop (a) and the 
distribution of current during the ascending (b) and descending branch (c) of the loop. The 
arrows indicate the increasing and decreasing of current, i.e. the time evolution. The dot 
lines represent the boundary between S and W regions. ds/dvi/=l/2, dl/dt=lA/s, jcTy=2x 
10«A/m2. 

(2.1) ]cs/'}cw =10, high inhomogeneity of flux pinning. Average speed of current diffusion 
in S region is too small comparing with that of W region, resulting that in S region there 
is no enough time for current diffusion in response to both the ascending sweeping and the 
descending one at the given parameters. No obvious V-I loop. 

(2.2) ]cs/']cw =4, medium inhomogeneity of flux pinning and lager average speed of 
current diffusion in S region, resulting that more current has penetrated into the S region 
during the descending branch than the ascending one at the given parameters. Obvious V-I 
loop is observed. 

(2.3) ']cs/']cw =1, uniform flux pinning, equal average speed of current diffusion in the 
two regions. There is little difference between the distributions of current during the two 
branches. No obvious V-1 loop. 

Fig. 3. Numerical results, showing the effect of dg/'^w ratio on the V-I loop (a) and the 
distribution of current during the ascending (b) and descending branch (c) of the loop. ]cS 
= 4 dl/dt=lA/s, jeH^=2xl08A/m2. 

(3.1) d^/dvK =0.857, the thickness of S region is proper, resulting that in S region there is 
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enough time for current diffusion in response to current descending, but the current diffusion 
is too busy to respond to current ascending for the given parameters, and area of the V-I 
loop is large. 

(3.2) dg/dvr =0.222, S region is thinner and the distance of current diffusion in S region 
is shorter, resulting that in S region there is much more time for current diffusion in response 
to current ascending at the given parameters, and the area of V-I loop becomes smaller. 

(3.3) ds/dw =0.105, S region is too thin and the distance of current diffusion in S region 
is too short, resulting that in S region there is enough time for current diffusion in response 
to both the current ascending and descending at the given parameters, no obvious V-I loop. 

Fig.4. Numerical results. Current sweeping rate dependent V-I loop (a). The corresponding 
distributions of current during the ascending (b) and descending branch (c) were also shown. 
ds/dw=l/2, ],s/]cw =4, jcH^=2xlOSA/m2. 

(4.1) dI/dt=O.OlA/s, the current sweeping rate is so slow that in S region there are 
enough time for current diffusion in response to both the ascending and descending current 
sweeping at the given parameters, no obvious V-I loop. 

(4.2) dI/dt=2A/s, the current sweeping is faster, resulting that in S region there is no 
enough time for current diffusion in response to current ascending at the given parameters, 
and the V-I loop is larger. 

(4.3) dI/dt=10A/s, the current sweeping is so fast that in S region there is no enough 
time for current diffusion in response to both the ascending and descending current sweeping 
at the given parameters. No obvious V-I loop. 

Fig. 5. Numerical results. Magnetic field dependent V-I loop, where the area of the loop is 
insensitive to }cW while the transport jc is strongly affected by magnetic field. d5/di^=l/2, 
jc5 =8xl0SA/m2, dI/dt=lA/s. 

Fig. 6. Experimental V-I loops of the Ag-Bi2223 tape with different current sweeping rates 
dl/dt, which is qualitatively with the numerical curves shown in Fig.4 and is a conformation 
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of our model. 



Fig. 7. Experimental V-I loops of the Ag-Bi2223 tape at different applied magnetic fields, 
which is qualitatively with the numerical curves shown in Fig. 5. It is noted that the area 
of the loop is insensitive to Ba whereas the value of the jc is strongly affected by magnetic 
field, which supports the present model further. 
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